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We present a study of the cross sections J/tp X — » D^*' D^*' (X = p, $) based on the calculation 
of the effective tri- and four-linear couplings J/ip(X)D^D^ within a constituent quark model. In 

ry ■ particular, the details of the calculation of the four-linear couplings J/tpXD^*' D^*' are given. The 

results obtained have been used in a recent analysis of J/tp absorption by the hot hadron gas formed 

s ! . in peripheral heavy-ion collisions at SPS energies. 
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CN I I. INTRODUCTION 

o : 

io ■ 

The problem of computing J/tp strong couplings to ir, p and other pseudo-scalar and vector particles has its own 
' interest because it opens the way to the calculation of cross sections of the kind (see Fig.^): 

'-H ; ^(j/vk/v..} — > dW5w). a) 

Such cross sections are the basic ingredients to estimate the hadronic absorption background of J/tp in heavy-ion 
collisions, as it is thoroughly discussed in 0,0- The description of processes like those in Eq. JTJ is a hard task 

. ' because no experimental test can be performed and moreover they are not amenable to first principles calculations, 
. so that one has to resort to build models and make approximations to describe their dynamics. 
;_i ' The dissociation process of the J/tp by hadrons has been considered in several approaches, but the predicted cross 

5^ i sections show very different energy dependence and magnitude near threshold. Anyway, using different approaches, 
one consistently finds non negligible cross section values (at least comparable with the nuclear one AfJ/tp, Af =nucleon) 
especially for the reactions with 7r's and p's, the most studied cases; for a review see for instance Ref. This is 
certainly a clear indication that the picture of J/tp absorption by nuclear matter, as an antagonist mechanism to the 
plasma suppression, is incomplete as long as interactions with the hadronic gas formed in nucleus- nucleus collisions 
are not considered. 

The problem of calculating the J/tp dissociation by pseudo-scalar and vector mesons has been addressed in Refs. 0,0 
within the Constituent- Quark-Meson model (CQM), originally devised to compute exclusive heavy-light meson decays 
and tested on a quite large number of such processes [j]. The basic calculations refer to ir and p contributions. The 
couplings to other mesons have been obtained under the hypothesis of flavour/octet symmetry. 

The typical effective Feynman diagrams contributing to the J/tp dissociation are depicted in Fig. ^ The tri-linear 
couplings pD^D^ have been calculated in and the J/tpD^D^ couplings have been recently discussed in (we 
report these results at the end of Sect. II), where also four-linear couplings involving pions have been derived. The aim 
of the present note is to explain the method used and the results obtained in evaluating the four-linear couplings of the 
kind J/tppD^D^ (see Fig. 1, third diagram), since these are not calculated elsewhere within the CQM framework. 

The numerical values of the J/tpQDs ' D\ ' couplings are also given. For completeness we report also the expressions 
for the tri-linear couplings discussed in Refs. 0,0- In the end, we present the cross section predictions, based on the 
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FIG. 1: Tree level effective Feynman diagrams for the J/ip p — > if if reaction, if being Z)'*', with = D or D*. 



complete set of contributing diagrams, for the processes J/ip p — > £)W flW and J/^ $ — > L'i*' 1 together with 

an estimate of the associated theoretical uncertainties. 



II. THE MODEL 



CQM is based on an effective Lagrangian which incorporates the heavy quark spin-flavor symmetries and the chiral 
symmetry in the light sector. In particular, it contains effective vertexes between an heavy meson and its constituent 
quarks (see the vertexes in the r.h.s. of Fig. EJ) whose emergence has been shown to occur when applying bosonization 
techniques to Nambu-Jona-Lasinio (NJL) interaction terms of heavy and light quark fields pj- On this basis we 
believe that CQM can be considered as a quite reasonable approach to the computation of J/ijj strong couplings to 
be compared to the various methods available in the literature, often based on St/ (4) symmetry 9j. 

In Fig.|2we show the typical diagrammatic equation to be solved in order to obtain 94(53), four(tri)-linear couplings, 
in the various cases at hand: on the l.h.s. it is represented the effective four-linear coupling to be used in the cross 
section calculation (to obtain one of the relevant tri- linear couplings we could discard either the J/ifi or the p); the 
effective interaction at the meson level (l.h.s.) is modeled as an interaction at the quark-meson level (r.h.s. of Fig.[5J). 

The J/ip is introduced using a Vector Meson Dominance (VMD) Ansatz: in the effective loop on the r.h.s. of Fig.0 
we have a vector current insertion on the heavy quark line c while on the l.h.s. the J ftp is assumed to dominate the 
tower of 1~, cc states mixing with the vector current (for more details see @). Similarly, vector particles coupled 
to the light quark component of the heavy mesons p, lj, when q — (u,d), or K*,$>, when one or both light quarks 
involved are q = s, are also taken into account using VMD arguments. 

The pion and other pseudo-scalar fields have a derivative coupling to the light quarks of the Georgi-Manohar 
kind Hg. 

In this paper we will mainly focus on the reaction: 

J /ip p — > D W f5 ( * } (2) 

and in particular on the four-linear coupling J/ippHH (third graph) in Fig. ^ 

In CQM, as in Heavy Quark Effective Theory (HQET) ^l]]> the heavy super-field H describes the charmed states 
D and D*, respectively associated to the annihilation operators P5, P M . H is written in the following way: 

H(v) = ±±l(r-P 5l5 ), (3) 

where v is the four-velocity of the heavy meson; the limit of very large heavy quark mass is understood. The heavy 
quark propagator is: 



l + i) i 
2 vk'' 



(4) 



where k is the residual momentum defined by the equation Pq = rriQV^ + k^ and related to the interaction of the 
light degrees of freedom with the heavy quark (k ~ O(Aqcd))- 

The p is described by the interpolating field p^ 7] and its kinetic term in the effective Lagrangian is built out by 
the tensor field strength: 

Fnv(p)=dnPv-dvPn + [Pn,P*>]- (5) 

In the approach followed by [7j, is defined by p^ = im p / f p p^ where p is a 3 x 3 hcrmitian traceless matrix analogous 
to the 3 x 3 7r matrix of the pseudo-scalar octet. 
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FIG. 2: Basic diagrammatic equation to compute the <?4 couplings. The l.h.s. is the effective vertex J/ippHH at meson level (e 
and r\ are respectively the p and J/tp polarizations); while the r.h.s. contains the Hoop process to be calculated in the CQM 
model. 



In the following we will use the Feynman rules defined in 4] . The interaction terms relevant to this calculation are: 




-qHQ v + h.c, (6) 
which describes the vertex light quark (q), heavy quark (Q v ), heavy meson (H), and 

(7) 

describing the vertex light quark, light quark, p. Here the decay constant f p is defined by 

(0\VMl,e))=ifpt» (8) 

where /„ = 0.152 GeV 2 . We use a similar definition for the decay constant of J/4> but we factor out the mass mj in 
the latter case, thus obtaining fj = 0.405 GeV. 

Once established the form of the effective vertexes occurring in the loop diagram in the r.h.s. of Fig. |3 one has 
just to compute it using some regularization; we will adopt the Schwinger proper time. 

CQM does not include any confining potential and an infrared cutoff p is needed to prevent low integration momenta 
to access the energy region where confinement should be at work. The kinematic condition for the free dissociation 
of H in vriQ and m q is given by: 

m H > m Q + m q (9) 

with ph = tuh v ~ toq v + k. It follows that 

k ■ v > m q . (10) 

In the hadron rest frame we have k° > m q and we can therefore require: 

pKm q . (11) 

The value of the constituent light quark mass in the model at hand is given by a gap equation £| : 

m g -m -8Gh(m 2 q ) = 0, (12) 

where G = 5.25 GeV -2 , mo is the current mass and the I\ integral is defined in the Appendix. As a consistency check, 
putting a zero current mass for the u, d species we get a constituent mass of 300 MeV, while for a strange current 
mass of mo = 131 MeV we obtain a strange constituent mass of 500 MeV using p = 300, 500 MeV respectively in the 
calculation of I\ . 

The residual momentum has an upper limit given by the chiral symmetry braking scale A x ~ 4^/^ which we adopt 
as a UV cutoff 0. 

Then the momenta running in the loop are limited by two cutoff's: p and A. These two cutoff's are implemented 
by the Schwinger regularization on the light propagator as follows: 

d 4 l n9 1 — > ( dH / lAl ds e~< l2+ <). (13) 
(P-m 2 ) J J 1/A 2 
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The diagrammatic equation in Fig. [3 states that the effective vertex J/ippHH is given by: 



(-1) \J Z H m H Z H >m H , 



dH 



v' ■ I + A /,/mj u • / + A y " 1 - m f p f + $-m 



(14) 



H and H' represent the heavy- light external meson fields labeled by their four- velocities v,v' while the 
V Zn'nriHZH'fTiH' coupling factor of heavy mesons to quarks is calculated in The parameter A appearing in 
the heavy propagator is defined by: 

A = Mji — rriQ, (15) 

i.e., the mass of the heavy-light meson minus the mass of the heavy quark contained in it. A is the main free 
parameter of the model. It varies in the range A = 0.3,0.4,0.5 GeV for u, d light quarks and 0.5,0.6,0.7 GeV for 
strange quarks Varying A gives an handle to estimate the theoretical error, m is the constituent mass of light 
quarks as defined above. 

III. THE CALCULATION 

The p is coupled to the light quarks by VMD, e being its polarization and q its 4-momentum. The J/tp, having 
polarization rj, is also coupled via VMD, but to the heavy quarks (rj appears in the trace between the two heavy 
quark propagators, while e appears between the two light quark propagators). In front of this expression we have the 
fermion loop factor. 

The trace computation in l|14|l will introduce a number of scalar combinations of the momenta and polarizations 
of the external particles that we will list in the Appendix. Each of these combinations will be weighted by a scalar 
integral which amounts to a numerical factor: what we call the coupling. Actually, as we will see, such scalar integrals 
depend on the energy of the p. In general the expressions obtained for the four-linear couplings appear to be quite 
complicated functions of E p if compared, e.g., to those obtained when studying only J /if; interactions with pions [||. 
It is therefore difficult to write down general polar expressions for the E p behaviour. On the other hand we have 
in mind to use these results to compute cross sections (Jj/^ p and thermal averages (p ■ <j)t in a hadron gas at a 
temperature T w 170 MeV where the Boltzmann factor is presumably more effective than any polar form factor in 
damping the high energy tails. 

Using the Feynman trick the fermion loop of the above equation i|14fl becomes: 

m J m P /y 7 [\ 9 -at f dH Tr [H' j H (1/ - jx + m) j (f - jx + j + m)\ 

7^7 T P VZHmH ZH,mH 'L dX d^) lNc J (/ 2 - fa 2 ) (vl + S) + ' (16) 

in which we have defined: 

m 2 — m 2 + x m p (x — 1), (17) 

and 

S = A- xq-v = A- xE p (18) 

6' = A- xq-v' = A-xujE p , (19) 

where E p is the energy of the incident p and lj — v ■ v' (v' = ljv) . The cross section computation is performed in the 
frame where J/ip is at rest, uj, in this frame, is related to the meson masses by: 

_ ™j/v- + m l - m H- m w + 1E p mj/^ 
lu — . (20) 

By kinematic considerations the energy threshold of the reactions (j2J) for DD and D*D channels is E p ~ 0.77 GeV 
whereas for D*D* channel is E p ~ 0.96 GeV, with wwl. We consider p particles with energies in the range between 
0.77 and 1 GeV where the two final state mesons are almost at rest. 

All the couplings that we can extract by direct computation can be written in terms of 7 basic expressions which 
we call: L 5 , A, B, C, D, E, F. The latter are linear combinations of the ij, L t integrals listed in the Appendix and are 
defined by: 

,f!n, ^w.. F 1 ,w.., , , ^ =L 5 (21) 



dm 2 J (2tt) 4 (I 2 -fa 2 ) {v I + S) («'•/ + 5') 
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dm 2 iNc J (2^) 4 {P -rh 2 )(v-l + S) (v'-l + 8') ~ A " M + B ^ (22) 
iN ° / 7^1772 - 2w rf^w^ n-^ = C 'g^ + J> ^+^<< + Ji1 (^X+<^)- (23) 



d?7i 2 7 (2tt) 4 (Z 2 - m 2 ) («•/ + <*) (v' • Z + <5') 
The final expression of the loop integral can the be reduced to a sum of terms of the general form 



J2S(H,H') C [ dx gf\x,E p 

a JO 



(24) 



where S{H 1 H t ) represent the scalar combinations of momenta and polarizations of H and H l occurring in the 

(S) 

calculation; g\ , are the corresponding couplings. Here C is given by: 

771^ 772 ^ 

C = -A y/Z H m H Z H ,m H ,, (25) 

Jjmj f p 

with fj = 0.405 GeV. Our central values for the couplings are obtained for A = 0.4 GeV (and Zh = 2.36 GeV -1 Q), 
while for thh we use the experimental value for the mass. 

In Table|Uwe list the explicit expressions of the couplings. We call them <?4 = {g, h, /}, respectively related to the 
four linear couplings J/ippDD, J/ifjpD* D, J/tppD*D* . All these expressions have to multiplied by C and integrated 
over x; the numerical values in Table I are then: 



C l dx 9i{x,E p ). (26) 
'0 



These are typically complicated functions of E p and it is not as easy as in [6j to determine an explicit polar form 
factor dependency common to all of them. On the other hand we have in mind to adopt these couplings to compute 
cross sections of the kind crj/^ p ^ p Cn charm and use this information to compute thermal averages in the form: 

\P' cr J/i/>p-+opcn charm/ T — ^2 dE g£ ;/ fcT _ J i (27) 

where Eq is the energy threshold needed to open the reaction channel and the Bose distribution is used to describe 
an ideal gas of mesons, p in the l.h.s of Eq. (|27|l is the number density of particles in the gas. The Boltzmann factor 
exp (-E/T) will be at work as an exponential form factor cutting high energy tails faster than any polar one. We 
could therefore avoid any arbitrary Ansatz on form factors at the interaction vertexes. We limit ourself to study the 
dependency of our couplings on E p in the range of energy where we reasonably think to have p mesons in the hadron 
gas excited by a peripheral heavy-ion collision. Estimating the J/ip absorption background to the suppression signal 
in heavy- ion collisions amounts to compute the attenuation lengths (inverse of the thermal averages in Eq. (|27|0 in a 
hot gas, T w 170 MeV, populated by n, K, T], p, to, ... mesons. We could therefore expect to have E p w 770-^ 1000 MeV. 

The loop that must be calculated in the case in which we substitute a $ particle to the p is the same as in Fig. [5] 
but with q,q' — s (m s — 500 MeV, p ~ fn s ), and with super-fields H s in place of H. The reaction in this case is 

J $ — ► 5W. (28) 

The structure of the coupling of <& to light quark current is identical to the one of the p, and all the above equations 
are valid also in this case with the substitutions: m p — > m$, f p ~*f^,H—f H s . Numerically we have used /$ = 0.249 

GeV 2 , while for m# s we have used the experimental value for Di*\ The numerical values are reported in Table [I] 

We conclude this section by writing the tri-linear couplings pHH; they are computable within the same framework 
by simply not including the J/ip in the diagrammatic equation of Fig. [21 

The vertex pD^D^ is described by two constants g% = f3, A and the effective Lagrangian describing this interaction 
can be written as 0, Q : 

Chh p = -ifiTrlHH] v ■ p + iX r Tr[Ha flu H]J r ^, (29) 
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J/ipXD*D* X = p X = $ 



/i -^-(A + B -mL 5 ) 25.5 ± 0.6 13.6 ± 0.1 GeV" 2 

/ 2 ^rr( B -mU + A(2lux - 2x + 1)) 26 ± 1 13.86 ± 0.08 GeV" 2 

f 3 dh Ax 4 ± 2 1.5 ± 0.5 GeV" 4 

fi {A - mL 5 + b1-2lux + 2x + 2lu - 1)) 26.0 ± 0.5 13.8 ± 0.1 GeV" 2 

/s ^-B{x-1) -2.2 ± 0.7 -1.2 ± 0.1 GeV" 4 

/e {-m 2 L 5 + 2C + D + E D +2Fuj + m 2 p (x 2 - x)L s )[l - u) 0.03 ± 0.1 0.03 ± 0.03 

f 7 7; ^-(~m 2 L 5 + 2C + D + E + 2Fuj + m 2 p (x 2 - x)L 5 ) -0.1 ± 2 -0.2 ± 0.2 GeV" 2 

fs ' -J-(D + F -mA) -8 ± 2 -8.2 ± 0.3 GeV" 2 

f 9 -^—(—m 2 L5 + 2mA + 2C — D + E + 2F(lj — 1) + m 2 (x 2 — x)L$) 8.3 ±2 8.0 ± 0.1 GeV" 2 

ho 7 ^-(-m 2 L 5 +2mB + 2C + D-E + 2F(Lu-l)+m 2 p (x 2 - x)L 5 ) 7.3 ± 0.8 6.3 ± 0.2 GeV" 2 

fn -4-{E + F-mB) -7 ± 1 -6.5 ± 0.5 



TABLE I: The couplings g± = { gi, hi, fi} expressed as linear combinations of the basic scalar integrals listed in the Appendix. 
The numerical values are given by C J Q dxgA '■ the mean values are estimated by setting A = 0.4 GeV (A = 0.6 GeV) and 
varying the energy of p ($) in the range E p — 0.770 — 1 GeV, (E<s> = 1.02 — 1.2 GeV), while the error is calculated by combining 
the excursion in the selected energy range at A = 0.4 GeV with the uncertainty over A (in the range A = 0.5 — 0.7 GeV). 
Some of the couplings are unavoidably affected by large uncertainties. This affects the determination of the attenuation lengths 
discussed in the text to the extent pointed out in p|- 



where the field p M and the tensor T^" have been defined above. The numerical values are: 



(3 = -0.98 (30) 
A = +0.42 GeV" 1 ; (31) 



7 



analogously /3 and A for the Chh® are 



(3 
A 



-0.48 

-0.14 GeV" 1 . 



As for the couplings J/ipD^D^*' , they have been extensively discussed in 
Observe that 

£j/4,hh = igj/i/,HHTt[H'yf l H]J li , 



(32) 
(33) 

|. Here we just report the main results. 



(34) 



where H can be any of the pairs D D* or D S D* (neglecting SU(3) breaking effects). As a consequence of the spin 
symmetry of the HQET we find: 



9j/ipD*D' — Sj/ijjDD i 



9j/i>DD* 



g.j/4>DD 



(35) 



The numerical values are given by: 



g.j/ipDD = 8.0 ±0.5 
gj/ifiDD* = 4.05 ±0.25 GeV- 1 
gj/i>D*D* = 8.0 ±0.5. 

In Fig. we report the cross section curves as functions of the y/s of the process for the three final states under 
consideration (DD, DD* , D* D*). This calculation has been made by using the tri- and four-linear couplings quoted 
above, assuming their validity in the energy range y/s « 3.8 -j- 4.5, and computing the tree level diagrams for the 
process at hand (for a sketch of the diagrams involved see pj). The dashed curves define the uncertainties bands 
obtained by varying A and E p , as discussed in Tab [I] 
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FIG. 3: The cross sections of the processes J/ip p — ► and J/tp $ — > Da D s *' on the left and on the right panel 

respectively. The dashed curves define the uncertainty bands obtained by varying A and E p as discussed above. Some of the 
reactions, the one initiated by p giving DD in the final state and those initiated by cj> giving D 3 D S and D*D S (or D 3 D*, a sum 
of the two is taken) show the typical "exothermic" peak for zero p((j>) momentum. The remaining reactions show the usual 
threshold behaviour (endothermic). 



IV. SUMMARY 



We have presented the calculation method of the effective couplings J/tp(X)D^ D", with X = p, within 
the CQM model. The resulting cross section predictions, together with an estimate of the associated theoretical 
uncertainties, have been presented as functions of y/s, showing values of the same order as the cross sections for 
J/V>tt -> £)(*)£)(*). This, given also the higher spin multiplicity of the p meson with respect to pions, demonstrates 
the importance of the p contribution to the J ftp absorption in the hot hadron gas, formed in peripheral heavy-ion 
collisions at SPS energy, as discussed thoroughly in Rcf. [2j . Aiming at calculating thermal averages with T w 170 MeV, 
we didn't discuss in the present paper the introduction of any arbitrary form factors since the exponential statistical 
weight acts as a cut off in the high energy tail. 



8 



Acknowledgments 

We wish to thank L. Maiani for the stimulating collaboration and encouragement. 



Appendix 

In this Appendix are listed the Ii and Li integrals occurring in the calculation and their linear combinations 
A,B,...,F. These integrals have been computed adopting the proper time Schwinger regularization prescription, with 
cut-off \x = 0.3 GeV (0.5 GeV when is present a strange quark), A = 1.25 GeV. In the following N c = 3. 



h = iN c 



dH 



1 



(2tt) 4 (I 2 - m 2 ) 



N, 



c m 2 T | -1 



h(S) 



16tt 2 

-iN c 
N c 

16^3/2 J 1/A2 

dH 



A 2 ' /j, 2 
dH 1 



(36) 



(2tt) 4 (I 2 - m 2 ) {vl + S) 
^ ds 



L^e-^-n ( 1 + Erf(^)) 



(37) 



I 5 (5,S',lj) = iN c J ^- 



1 



dy-_ 



Y (l 2 - m 2 ) (vl + S) [v'-l + 5') 
1 



o l + 2y 2 (l-Lu) + 2y(LU-l) 

6 n/*" 

16^3/2 



+ 



/ " ds a e - s (* 2 - 2 ) s -i/2 (i + Erf(<7>/i)) 

Jl/A 2 

4 



ds e- sa s- 1 



in the last expression we have defined 

a = a(S,5',y,w) = 

In the previous equations m 2 , 6 and 8' are given by 



6(1 -y) +6' y 



y/l + 2(u-l)y + 2(1 -uj)y 2 ' 



~ 2 2 , 2 / i \ 

m = m + .x m p (x — 1) 

5 = A — x g • v = A — xE p 

S' = A — x q ■ v lli = A — xuj E p , 



(38) 



(39) 



(40) 
(41) 
(42) 



with m = 0.3 GeV the constituent mass for light quark u, d. The expression of uo = v ■ v' in the rest frame of J/ip is 



TO JM + m 2 p -m 2 H- m H> + tEp-mj/^ 



2mHmH> 



In the I\ integral the gamma- function is 



T(a,x 0} x 1 ) = dt t 

J x Q 



t j.a-1 



while the error function is 



Erf (a) = jL(dx e- x \ 
V 71 " Jo 



(43) 



(44) 



(45) 
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The Li integrals are denned in the following way: 



and they are 



U = 



dm 2 



u = 



L 3 (5) 



d 
dm 2 

N c 
16tt 2 

d 



iN, 



dH 



1 



(2tt) 4 {I 2 - TO 2 ) 

r ! -l 



~ 2 ~ 2 
TO TO 



-iN r 



A 2 ' 
dH 



fh 2 



- 2 ~ 2 \ ~\ 
TO TO x 



r -l, 

dm 2 \ ' A 2 ' /i 2 



1 



dm 2 J (2tt) 4 (Z 2 -to 2 ) (u ■/ + <*) 



16tt 3 / 2 



l/A 2 



e - s (™ 2 - 52 ) (-s- 1/2 ) (1 + Erf (5Vi)) 



_9_ 

dm 2 
6 



d 4 / 



1 



(2tt) 4 (/ 2 - to 2 ) (u • I + 8) (v'-l + 8') 



dy- 



1 



16tt 3 / 2 7 y l + 2?/ 2 (l-cj) + 2y(tj- 1) 
x d Sff r s ("' 2V ) f-s 1 / 2 ) (l + Erffcr^)) 

Jl/A 2 V 7 

The A, B, C... functions are all function of x and E p through 8, 8' and uj: 

L 3 {5') + 6> L 5 (8 7 6', w) - (L 3 (5) + <J' L 5 (S, 8', u)) u 



A 
B 



uj 2 -\ 

L 3 (8) + 8' L 5 (S, 8', uj) - (L 3 (8') + S L 5 (S, 8', uj)) uj 
uj 2 - 1 



(46) 



(47) 



(48) 



(49) 



C 



D 



1 



2(uj 2 - 1) 



L 5 (<5, w) <5' 2 + (i 3 (*) " (i3(0 + 25L 5 (<5, <5', «)) + J(L 3 (<J') + ££ 5 (<5, 5', «)) 



-<5 L 3 (<5) w + I 5 (<5, w) - 1) + ^ 5 (<5, 5', w) ™ 2 (^ 2 - 1) 

' 2(Li + 5i 3 (,5)V 3 + (4(<5, 5', w) + 2<5'(L 3 (5) + <5'L 5 (<5, 5', w)) + L 5 (8, 8', uj) m 2 )uj 2 ~ 2L x uj - 5SL 3 (8)uj 



2(uj 2 - l) 2 



-I 5 (S, 8', uj) + 8' 2 L 5 + 3S(L 3 (S r ) + 8 L 5 (5, 8', uj)) - L 5 (S, 8', uj) to 2 + S'(L 3 {6) - 3{L 3 (S') + 2SL 5 (5, 8', uj))uj) 



E = 



1 



2(Li + 8'L 3 {8') )uj 3 + (I 5 (6, S', uj) + 28{L 3 {8') + SL 5 (8 7 8' , uj)) + L 5 (S, 8', uj) m 2 )uj 2 - {2L X + b8'L 3 {8') 



2{uj 2 -\) 2 

+3S(L 3 (8) + 2S'L 5 (8, 8', uj)))uj ~ h{8, 8', uj) + 3S'L 3 (S) + SL 3 (S') + 8 2 L 5 (S, 8', uj) + 3S' 2 L 5 (S, 8', uj) - L 5 (S, 8', uj)m 2 
F = of 2 1 ^2 f(-( J 5(^ w) + L 5 (S, 8', uj) to 2 ) V 3 + {-2L X + 8'L 3 {8') + 8(L 3 (8) + 4S'L 5 (S, 8', uj)))uj 2 

2{uj a — \y [ 

+ (I 5 (S, 8', uj) - 3(8(L 3 (8) + 8L 5 (8, 8', uj)) + 8'(L 3 (8) + S'L 5 (8, S\ uj))) + L 5 (S, 8', uj) m 2 )uj 
+2L 1 + 2(8'L 3 (S) + 8(L 3 (8) + 8'L 5 (S,S',uj))) . 

The determined couplings weight the scalar combinations of external particle momenta and polarizations in the 
combinations obtained by the loop computation in Eq. (|14|) . We list below the scalar combinations S 
together with the couplings they are weighted by, and with the effective Lagrangian interactions £ one can write 
down for them. These terms describe the four- linear diagrams in Fig. 1; for the actual cross section computation one 
has to add also the t and u channels. 
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S(DD) 


n 

y 


£ 


Q 






-J a d^Pv d»D d v d a D 


q 


■ P2 V ■ P2 £ ■ Pl 


92 


J a d^Pv d^d a D d v D 


q 


■ pi e ■ p 2 V ■ Pi 


9i 


-J a d^d a D d u D 


Q 


■ pi e ■ p 2 rj ■ p 2 


-.92 


J a d,jp v d^D d u d a D 




rj-qe-pi 


53 


-d^ Pv J» d v D D 




VPi e 'Pi 


9i 


-J^p v d»d v D D 




VP2 £ ■ Pi 


95 


J-dD p-dD 




q ■ r] e ■ p 2 


-9e 


d^Pu J» d v D D 




q-Pi e-r) 


-93 


-d^ Pv r d»D D 




q ■ P2 e ■ 77 


96 


d^ Pu ,r d»D D 




e • P2 V ■ pi 


97 


J-dD p-dD 




C-P2VP2 


9s 


-j^p v d»d v D b 




V ■ e 


99 


J-pD D 



S(DD*) h C 



iSa/3-fS q a e l3 rf 


4 




hi 


-£a0j5 d a pl } J^D* S D 


i£a/3 7 S q a ^rfip{ 


m 


■P2 


1 h 2 

m D *m D 


-£ a p-fS d a pPj^d 5 D^d^D 


ieafSjS q a e ti tfpi 


m 


■P2 


\-h 2 


£ aPl 5 d a p^JW s d ll DD* 1 " 


iSap-yS q a e fi rjlp[ 


v- 


Pi 




£ a ^s d a pPd 5 d»D*~<J^D 


iSap-yS q a ePrftp\ 


rj- 


P2 


1 h 3 

m D *m D 


£ a01 s d a p^d 5 D* 1 J ■ 3D 


ie a f3~(6 q a e l3 i]Jp2 


v- 


Pi 


1 h A 

m D *mo 


-£ a ^8 d a p^d^D* 7 d s DJ^ 


ie a py8 q a e f3 r]Jp s 2 


v- 


P2 




£ a01 s d a pPD* J d s d^DJ» 


i£ af 3j8 q a e l3 pjp s 2 


v- 


m 


1 h 5 


-ZafS-yS d a pPd~>D^d S DJ^ 


iSap-yS q a rf > r]\p\ 


e ■ 


Pi 


(/i 5 - /i 6 ) 


£ a 0-yS d a p^Pd 5 d^D*~<D 


iSa/3-fS q a r) ViPi 


e ■ 


P2 




-£ af3l 8 d a p^Jf'd 5 D* 1 d ll D 


ie a py8 q a V f3 ViP2 


e ■ 


Pi 


1 K 


-£ a 0js d a p^J^d^D* 1 d s D 


iSa/3-fS q a V^ViP 5 2 


e ■ 


P2 




£a/3 7 5 d a p»jPD*~<d S d^D 


i£al3jS q a V f3 pJP2 


e ■ 


Vl 


1 h 7 

m D *mo ' 


-£ a lS d a p^J^d^D* tl d & D 


ie a [3~f5 ^rfrilPx 


Q- 


Pi 




£ a(il 8 d^ P a Ji 3 d s d fi b* 1 D 


i£af3-fS ^rfnlpi 


q ■ 


P2 




-£af3 7 5 d»p a jPd 5 D*~'d tl D 


i£ a /3 7 s e a 'n l3 ViPi 




r h 2 

m D *m D 


-£ aPl s p a J d s D* 1 D 


ie a f3~(6 e a rPr){p\ 




Pi 


-£ a(il 8 d^p a Jf ) d^D* 1 d s D 


i£ap-fS e a V f3 riiP2 


Q- 


P2 


^h 2 

m D 


£ a(il 8 d"p a J^D* 1 d 6 d li D 






— -h 8 


-£ aPl 8 p a J fi b* 1 d s D 


ieaPjS e a ri l5 p'lp s 2 


q- 


m 




-£ af 3^8 d»p a JPd^D%d 5 D 


a 13 7 5 

i£ a fi~f8 e V1P1P2 


q 


■ V 


1 (h 2 + h 6 ) 


-£ af)l 8 d^p a d^D* f3 d s DJ^ 


a 7 8 

i£ a p 1 8 V V1P1P2 


e ■ 


pi 


— 3 — " hg 

m^m D 


-£ aPl s J a d^d^D* P d s D P ^ 


i£ a p 7 8 ^ViPiPl 


e ■ 


P2 


~ -i hi n 


£ aPl 8 J a d^D* P d 5 d^Dp^ 
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